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Abstract 



The transfer matrix method is applied to finite quasi-lD disordered samples attached 
to perfect leads. The model is described by structured band matrices with random 
and regular entries. We investigate numerically the level spacing distribution for 
finite-length Lyapunov exponents as well as the conductance and its fluctuations for 
different channel numbers and sample sizes. A comparison is made with theoretical 
predictions and with numerical results recently obtained with the scattering matrix 
approach. The role of the coupling and finite size effects is also discussed. 
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I. INTRODUCTION 



Since the discovery of non-self-averaging conductance fluctuations in mesoscopic conduc- 
tors Jl|-0, it has become necessary to reconsider various fundamental aspects of quantum 
electronic transport. In particular, the one-parameter scaling hypothesis fij has been chal- 
lenged || on the ground that the conductance G is not a self-averaging quantity, and thus 
it is necessary to consider its entire distribution rather than just the mean value. One of the 
approaches that provide the possibility of obtaining the conductance distribution is based 
on the theory of random matrices (for a recent review see, for example, ||). This theory 
exploits the connection between the conductance and finite-length Lyapunov exponents ]7[ 
of the transmission matrices (see also || and references therein) . Thus, in this approach, the 
main interest is in the properties of the Lyapunov spectra (all Lyapunov exponents ordered 
in an increasing/decreasing way). 

In this paper we study quasi-lD disordered samples of finite size described by Band 
Random Matrices (BRM). Such matrices have been extensively investigated in applications 
to thin wires (see the review || and references therein) or, equivalently, to ID Anderson-type 
systems with long-range random hopping. Scaling properties of the localization length of 
eigenfunctions are already well understood since some years |T0|j9[1 . More recently, the scaling 
properties of averaged Lyapunov exponents have been also established both for infinite [TT 



and finite [12| BRMs. In this paper, we are mainly interested in statistical properties of the 
Lyapunov exponents and thereby of the conductance as a function of the sample size and 
the number of transmission channels. 

One of the results of this paper concerns the comparison with the predictions of the 
Random Matrix Theory (RMT). It is known that, in contrast to other characteristics like 
the spectral rigidity, the shape of the level spacing distribution P(s) is quite insensitive to 
deviations from the Wigner-Dyson surmise, i.e. the prediction of the RMT. Nevertheless, in 
spite of an overall good correspondence between the P{s) numerically determined and the 
theoretical expression, here we show that the normalized second moment of the distribution 
P{s) reveals clear deviation from the theoretical predictions when applying the RMT to 
the Lyapunov Spectra. Accordingly, this deviation can serve as a sensitive measure of 
the correspondence of the data to the RMT results. We attribute such a discrepancy to the 
existence of a non-random repulsion acting on the low Lyapunov exponents from the left side, 
due to the symmetry of the Lyapunov spectrum. Though the observed deviation is relatively 
small, it manifests itself in the value of the variance for the conductance fluctuations which 
we have studied numerically by making use of the Landauer formula expressed in terms of 
the Lyapunov spectra. 

The outline of the paper is as follows. In the next Sec. II we describe the physical setup 
which permits us to study the properties of conductance of finite samples attached to the 
leads. First, we introduce the model in terms of Hamiltonian band matrices and derive the 
corresponding transfer matrices. Using these matrices, one can compute the conductance 
of quasi- ID finite samples through the Landauer formula. This is done by introducing a 
specific matrix, the eigenvalues of which represent the key ingredient of the conductance 
expression. At the beginning of Sec. Ill, we give a very short outlook of theoretical results 
about the connection between fluctuations of the conductance and of finite-length Lyapunov 
exponents. According to these findings, the statistical properties of Lyapunov exponents may 
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be approximately described by the RMT. In particular, the "level" spacing distribution is 
expected to follow the famous Wigner-Dyson expression. In this Section we compare these 
predictions with our numerical data obtained for the whole Lyapunov spectrum, finding 
clear deviations not only in the upper but also in the lower part of the spectrum, where a 
good correspondence with RMT was expected instead. A study of both the conductance 
and its fluctuations in different regimes, ranging from the very localized to the metallic one, 
is performed in Sec. IV. In particular, we have determined the dependence of the average 
conductance on the related scaling parameter when passing from the localized to metallic 
regime. Another quantity of interest was the universal conductance fluctuations for which 
we have checked the theoretical predictions. Detailed discussion of our results is presented 
in the last Section. 



II. THE PHYSICAL SETUP 

The simplest general model in the class of quasi-lD or ID systems with long range 
hopping is represented by the Schrodinger equation with interactions described by Band 
Random Matrices. This ensemble is defined as a set of real symmetric matrices the entries 
of which are independent Gaussian variables with zero average and variance ctq = 1 + 5 n ^ m 
(5 n>m is the Kronecker symbol) if \n— m\ < b, and zero otherwise (i.e., there are 26+1 nonzero 
elements in any row). The parameter b defines the hopping range between neighboring sites; 
in the quasi-lD interpretation, it is the number of transverse channels along a thin wire. 
In the study of the electronic conduction in disordered wires at zero temperature, one uses 
mainly one of the two theoretical approaches based on ideas of Thouless and Landauer 
respectively. In the Thouless approach, the wire is regarded as a closed system, and its 
conductance is defined through the sensitivity of eigenvalues to changes of the boundary 
conditions [p^-|15|. In the Landauer approach, one has to embed a disordered sample into 



a perfectly ordered lattice, i.e. to add perfect leads both to the left and to the right of the 
disordered sample. Thus, in contrast with the former approach, the latter one deals with an 



open system (for a comparison of the Thouless and Landauer conductance see, e.g., |L6| ) |l7|| ). 
Our model is described by the time-dependent Schrodinger equation, 

dcJt) n+b 

i 3^ = H n ^ m c m , (1) 

m=n—b 

where c n (t) is the probability amplitude for an electron to be at site n and H nm is a 
symmetric banded matrix with random entries for the disordered part (see below). The 
eigenvalues can be obtained by substituting the relation c n (t) = exp(—iEt)ilj n in Eq. (P and 
solving the resulting equation for ip n , 

I I n+6-l 

tpn+b = TT Eljj n - Yl H n,m1pm ] ■ (2) 

tln,n+b y m=n-b 

By defining x n {i) = ipn+b-i, the above equation can be recasted in the form of a 2b- 
dimensional linear map T n , 
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1 / 26 

X n+ i(l) = — I Ex n (b) -Y, H n,n+b-jXn{j) 



n,n+b \ 3=1 

= ar„0' + 1) K J < 26 (3) 

which provides an alternative interpretation of the stationary Schrodinger equation. In this 
picture, an eigenstate of Eq. ([I]) can be treated as a "trajectory" evolving under the action of 
map (S), and its localization properties are determined by the Lyapunov exponents. In the 



previous paper |T]J , we have investigated the shape of the Lyapunov spectrum in the limit of 
infinitely extended disordered samples. Here, aiming at a more complete understanding of 
conductance properties, we choose to rely on the multichannel Landauer formula 0] for the 
conductance of finite samples. The physical set-up requires considering a disordered sample 
of length L in between two perfect leads. At variance with the standard Anderson model, 
where only nearest-neighbor couplings are included, here the long-range hopping terms allow 
some freedom in the structure of the ordered leads and, especially, in the connection of the 
leads with the disordered sample. As for the leads, a natural way is to assume a band 
structure in the ordered part (leads) with H n)Tn = V, for \n\, \m\ > L (for the sake of 
simplicity we put V — 1). As for coupling matrix elements, we choose them in the same 
way as in the bulk, (random Gaussian entries with zero mean and variance cr 2 = 1). 

There is an immediate analogy with the Anderson problem on a stripe of width b, 
where one deals with lb x lb matrices, too. However, in our case the one-step ma- 
trix T n defining map (|3|) is not symplectic and the determinant is not equal to 1, but 
to (— 1) 26+1 (— Hb+i,i/ Hb+i,2b+i) ■ Yet, the total transfer matrix Tn = IT^Li^n, where 
N = L + 2b, satisfies the following relation 

(T N yZT N = £ (4) 

where T N denotes the Hermitian conjugate of Tn and the matrix E is defined as 

(5) 

Here S is a lower triangular matrix such that = 1, for any 1 < j < i < b and S* denotes 
the transposed operator. In fact, Eq. @ means that Tn has a "generalized symplectic" 
structure |18|| . 

The scattering properties of the sample can be better investigated by choosing an ap- 
propriate base, namely by decomposing the eigenfunctions in the plane waves, 

e inp 

= -7= , (6) 

V 2n 

supported by the ordered lattice. The corresponding eigenvalues are 

s i n f 26+i„) 

E = l + 2cosp + ... + 2cos(6p) = — . 2 , , (7) 

sm(f) 

while the velocity of a wave packet centered around p is 
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dE 2& + l cosggp) aiffl 

«(P) = ~T = — ^ ryh . ' , cos(p/2) . 8 

dp 2 sin(|) 2sm 2 (|) 

For simplicity, in what follows, we study the case E = for which all transmission channels 
are open. Therefore, there are b pairs of open channels with opposite velocities, 

?Trk 

P* = (-l)*26 + T > 1<I*I<*- (9) 

The matrix M connecting the amplitudes to the left with those to the right of the scatterer 
is 

' AR ^ = SjffifSo ( ) = M ( ) (10) 

where the matrix S is the transformation from the real space to the momentum space, i.e. 
it provides the expansion in momentum components. 

In order to determine the conductance, one needs to express the flux arising along the 
lattice. This requires a further similarity transformation 

F = TMT- 1 , (11) 

where the diagonal 2b x 2b matrix F is defined as T^j = Sij^/ui for i < b and Tij = Sijy/uiZb 
for i > b. The above transformation is equivalent to the normalization of the scattering 
matrix and it takes into account the fact that each open channel propagates with a different 
velocity. From Eq. (|H]) it is easily seen that the transformation (|TT|) corresponds to the 
change of variables af ,R = ^Jv^p^Af' 11 and b?"' R = y / t>(pj)B?"' R . The square of the modulus 
of the new variables gives the flux of particles which enter and exit the disordered sample. 
It is easy now to verify that F conserves fluxes, 

FV 3 F = a 3 (12) 

where 03 is a generalized Pauli a z matrix, i.e. 

and 1 denotes a b x b identity matrix. Condition ([12]) is crucial, since it corresponds to 
unitarity of scattering matrix. 

By using the symplecticity of F and F^F, and the relations between F and the scattering 
matrix, we find that the dimensionless conductance (measured in e 2 /h units) is readily 
written as M 



G = Tr 



p\F + (FtF)- 1 + 2 1 ^ 1 + cosh(2AT 7i (6, N)) 



S l + coshf2AM&,A0) ^ 



where the 7; (6, iV)'s are the Lyapunov exponents of the matrix F^F. 

The lengths 7 4 _1 (&, N) depend both on and on the realization of the disorder, but the 
ergodic multiplicative theorem Jl9| ensures that, in the limit N — > 00, i ~f^ 1 (b,N) converges 



towards the "localization length" l\ 



5 



lim j-\b,N) = l°°, 

V— >oo 



(15) 



i.e. lf° is a self-averaging quantity which depends on the energy E and on the bandwidth b 
only. Moreover, still in the limit N — > oo, the contribution of the similarity transformations 
introduced to pass from the transfer matrix T/v to F becomes negligible, implying that 
7i(6, N) are equal to the logarithms of the eigenvalues of the matrix 

O = lim [(Ttf)t(T*)] W*»> , (16) 

-L— >oo 

i.e. to the Lyapunov exponents of the product of the bare transfer matrices. 

In the following we shall denote with ji(b,N) > 72(6, N) > ... > 7* (6, iV) > ... > 
72b (&, iV) the effective ("finite-length") Lyapunov exponents |20] computed over a number 
L of iterations at energy E = 0. Because of the symplectic structure of F, the Lyapunov 
exponents ji(b,N) come in pairs with opposite values: for this reason, it is sufficient to 
compute only the positive exponents. 



III. LYAPUNOV EXPONENTS 

The standard Landauer approach to quantum transport in a two-probe geometry allows 
us to express the conductance G in terms of b real positive Lyapunov exponents 7$, see 
Eq. (pf). Thus, the knowledge of the joint probability distribution P({7i}) provides a 
complete statistical description of the conductance. Accordingly, we start this section by 
discussing the fluctuations of the Lyapunov spectra. 

In the metallic regime, an expression for P({7i}) has been proposed on the basis of 
Random Matrix Theory [pl-|2"4|. Namely, it was conjectured that the matrix 



jr=f tF + (^F)--2 (17) 

has the typical structure of full random matrices. This conjecture is based on the observation 
that F results from the product of many independent random matrices so that its entries 
are random numbers all of the same order. According to the RMT, the joint probability 
distribution of the Lyapunov exponents 7^ is of the type 

72, • • • , 7b) = C/9 exp[-/3#jv(7i, 72, • • • , 7b)] (18) 

where Cp is a normalization constant, (3 — 1, 2 or 4 is a symmetry parameter and is the 
effective Hamiltonian 

b b 

H N ( 7l , 72, • ■ • , 7*) = - E u ^ + E V(Tfi)- (!9) 

i<j i=l 

Eq. flT8D has the form of a Gibbs distribution at temperature T = for a fictitious system 
of classical particles on the semi straight-line (0, 00) under the external potential V(7i). As 
for the two-particle interaction, in conventional RMT it has the form 

u( 7ij 7i ) = - In I cosh(2A 7i ) - cosh(2A 7i )|. (20) 
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which corresponds to the logariphmic repulsion between the eigenvalues A« = cosh{2N r )i). 
All microscopic parameters are contained in the function V( 7 ), while the interaction poten- 
tial u is independent of them and has a geometrical origin (see the discussion in ||). The 
first impression was that, if ^(7) is suitably chosen, the distribution ( |l"8l , |T9"l) with (|20|) pro- 



vides an accurate description of transport properties in the metallic regime. However, recent 
developments |24| have shown that the expression (^0|) does not give the right value for the 
universal conductance fluctuations in the quasi-lD limit. The correct result is obtained if 
the two-body interaction potential (pO) is replaced by |25j (see details in 0), 



u{ lh 7j -) = ~ In I cosh(2iV 7l ) - cosh(2iV 7j )| - \ In | 7j 2 - 7 J| - ln(2JV) . (21) 

Although the above expression has been derived for the case j3 = 2 only, i.e. for broken 
time-reversal symmetry, there are indications that it is also valid for f3 — 1, 4. 

Less is known about the properties of Lyapunov exponents in the localized regime and 
in the metal-insulator transition (see || for a brief review). 

Below, we study numerically fluctuation properties of the Lyapunov exponents 7 i. Stan- 
dard techniques for the direct computation of eigenvalues are affected by large inaccuracies 
due to the high probability of small denominators in Eq. (^) . In order to avoid this problem, 
we have applied the algorithm originally developed in p6| for the infinite N— limit. The 
method consists in considering a formally infinite sequence of matrices all equal to F^F, 
which can be seen as "replicas" of the same disordered sample. The Lyapunov exponents of 
the corresponding product of matrices are then computed in a standard way by recursively 
applying the single matrices to b independent vectors which are continuously orthonormal- 
ized (see also p7|JT2]] ). The advantage of this approach is that one can keep the accuracy 
under control by renormalizing the vectors also in the intermediate steps that correspond to 
the application of the various matrices composing F. 

In practice the number of "replicas" has been fixed by imposing that the accuracy on the 
Lyapunov exponents is better than 10~ 4 for each disorder realization. Moreover, ensemble 
averages have been performed over sets of more than 2500 realizations. 

In the metallic regime, the RMT predicts that level spacings for eigenvalues Aj, and, 
therefore, the normalized differences 

1 < 1 < b ; "/,; (22) 



< 7i+1 - 7i > 



of the Lyapunov exponents are distributed according to the Wigner-Dyson (WD) surmise 
22| , |2"3|| . This expectation has been numerically confirmed in the two- [^2j and three- 



dimensional [27] Anderson model. 

In order to characterize the repulsion between any two consecutive Lyapunov exponents, 



it is useful to introduce the normalized width y of the "level" spacing distribution [27], 



( 2 3, 

where a denotes the corresponding standard deviation. For the WD-distribution, the pa- 
rameter y is equal to ywD = y 4/tt — 1 ~ 0.523. Values of y < ywD are interpreted as the 
signature of a stronger rigidity in the Lyapunov spectra. In particular, the limit y — > 
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corresponds to a (^-distribution, i.e. to a "crystallization" of the spectrum 0. For small 
values of y, the distribution of the spacings 5i around their mean values is quite close to 
Gaussian; this situation occurs in the localized regime when the repulsion is very strong (see 
below). Notice that in the complete absence of repulsion (which corresponds to a Poisson 
distribution), y is equal to 1 which thus represents an upper bound. 

It is known that the degree of localization in finite samples can be effectively described 
by means of the scaling parameter 

A = b 2 /N (24) 

which is connected to the b 2 proportionality of the localization length in infinite samples 
(see, e.g., [[| [K|). Numerical results for three different values of the scaling parameter 
A = 10.29, 6.95, 6.27 in the metallic regime are presented in Fig. 1 by open circles, diamonds 
and stars, respectively. The repulsion parameter y is reported versus the quantity \ = i/b, 
the standard parameter used in the representation of the whole Lyapunov spectra (see 
also [|TT],|TIJ). As one can see, the value of y is almost constant for x larger than 0.2. 
Although, in average, it is very close to the theoretical prediction ysw, there is a systematic 
difference, namely, y > ysw- This is a clear indication of a disagreement with respect to 
RMT expectations. The only exception is a narrow region around x — 1 where y is slightly 
smaller and thus more in agreement with the RMT. 

Larger Lyapunov exponents, in the region x < 0.2, are characterized by a much larger 
deviation from the WD-value. In particular, for x — > 0, the value of y decreases, revealing 
a much stronger repulsion than predicted by the Random Matrix Theory. This result is in 
contrast to previous studies for the 2D Anderson model p7| , |22| where it was found that the 



distribution for Si follows the WD-form for all Lyapunov exponents. 

Now, we analyze the statistical properties of the lowest Lyapunov exponent 7&(6, N) for 
different values of the scaling parameter A. Numerical data shows that the distribution of 
7h changes dramatically when A decreases. For example, for A m 0.6, the distribution V(Sb) 
turns out to be very different from the WD-form, and for A m 0.2 it looks like the Gaussian 
distribution around the mean value of 6b- In order to describe this transition quantitatively, 
we have fitted numerical data to the expression 



V(S) = A6?(l + B/36) f W exp -^(36 2 - £(1 - 0/2)6) (25) 




where f{(3) is the fitting function, f((3) = 2^(1 -0/2)/ (5- 0.16874. The constants A and B 
are defined through the normalization conditions 

J V{s) ds = 1; J s V{s) ds = 1 (26) 

while the free parameter {3 is the effective repulsion taking values in [0, oo). The above 
phenomenological expression has been suggested in Ref. |2S] (see for details) in order 
to describe the level spacing distribution for quasi-lD disordered samples of finite size, in 
the dependence on the degree of localization. A similar distribution has been also used in 
the description of the spacing between real particles (on a ring) interacting with each other 
via 2D Coulomb forces. In that case, f3 plays the role of the inverse temperature in the 
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thermodynamic equilibrium [^,^T[|. For (3 = 0, the dependence V(S) reduces to the Poisson 



distribution, while for (3 = oo it converges to a 5-function centered at 8 = 1. Moreover, for 
(3 = 1,2,4 the expression ( p5| ) is very close to that one given by the RMT (it even turns out 
to be more accurate than the WD-surmise, see ||29|| ). 



The "inverse temperature" (3 is a very appropriate indicator to characterize the transition 
from a totally uncorrelated Lyapunov spectrum (/5 = 0) to a "perfect crystal" {(3 = oo). 
In fact, it can be effectively determined by fitting the data with expression In order 

to emphasize the validity of our results, let us notice that the quality of the fit has been 
quantitatively checked by performing a x 2 test, that has been always brilliantly passed. In 
the metallic regime (e.g., for A = 10.24), the best fit of the numerical data for 75 gives 
(3 = 1 with a high accuracy, indicating that the minimum Lyapunov exponent is, indeed, 



characterized by the WD-distribution |23|. When moving towards the localized regime, 
where A < 1, the value of the repulsion parameter increases; for A « 0.2, we have (3 ~ 3.95. 
One can study the transition from one to the other regime by computing the repulsion 
parameter f3 for different values of b and N. The good data collapse observed in Fig. 2 
strongly indicate that (3 (for the minimum Lyapunov exponent) is a function of the scaling 
parameter A = b 2 /N only, a result which is in the spirit of previous studies fTOfl. Moreover, 
one can observe that already for A > 4, the repulsion parameter reaches the lower limit (3=1 
(we remind that, being the condition b N fulfilled, we are, for A ^> 1, in the metallic 
regime). For A < 1, the repulsion parameter starts increasing very fast and in the limit A — > 
it diverges. In other words, in the extremely localized regime, the probability distribution 
for the lowest Lyapunov exponent is a 5-like function, which can be approximated by a 
Gaussian. A similar effect has been observed in previous studies of 2D Anderson-type 
models [p7| , j32]1 . 

Finally, for what concerns the properties of the other Lyapunov exponents in the localized 
regime, we can look at full circles, diamonds and triangles in Fig. 1. First, we can notice that 
for all z's, including i = b, the quantity y(Si) decreases as we move towards the localized 
regime, leading to a Gaussian-like distribution. This is in agreement with the previous 
results for the minimum Lyapunov exponent 7& for which the same conclusion was reached 
by measuring directly the repulsion parameter (3. Moreover, in comparison with the metallic 
regime, the convergence to the 5-like distribution for the maximal part of the spectra becomes 
smoother. One should note that in the limit N — > 00, the distribution P(5i) for i < b was 
instead found to be Poissonian in the localized regime of the 2D Anderson model |27| , |3~2"ll . 
This means that the quasi-lD limit which we consider here, is somewhat special. Namely, 
the eigenvalues of the transfer matrix F*F are well separated, their number b being fixed 
while their mean separation increases linearly with N. On the other hand, if we consider a 
two-dimensional square of size N x N, then b will also diverge with N and the assumption 
of well separated eigenvalues is no longer valid. 



IV. CONDUCTANCE 

The present understanding of Anderson localization theory is mainly based on the one- 
parameter scaling. However, recent experimental results indicate that in order to understand 
the transport properties of disordered systems, one should study not only the mean expec- 
tation values but also the whole probability distribution. 
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One of the quantities that attracted a lot of research interests during last two decades, 
is the conductance. Theoretical and experimental results have shown (see references in ||) 
that the variance of the conductance is 0(1) (in e 2 /h units), independently of the average 
conductance, if the conductor is phase-coherent. A number of theoretical methods has been 
developed to compute the (universal) conductance fluctuations arising in the above men- 
tioned (metallic) regime. Many numerical calculations have also been performed although, 
until recently, mostly in connection with the Thouless approach which assumes that the 



system is closed (on this subject, see ]T7|]). Cm the other hand, very little efforts have been 
made to study open systems, taking into account the specific form of the coupling with the 
leads |16|,|18j, as in our case. Yet, it is known that the conductance of a mesoscopic sample 
depends sensitively on the measurement geometry and on the way the leads are attached 
to the sample |55|| . In this section we study the properties of the conductance G in the 
quasi- ID model ([!]) and compare them with theoretical expectations and numerical results 
obtained in a slightly different model [IK . 



Our first interest is in the geometric average of the conductance (InG) in the metallic 
regime 1 <C b <C N (here, the brackets (...) denote an average over different realizations of 
the random potential). As expected, the average conductance G = exp ((lnG)) turns out to 
depend on b and N only through the scaling parameter A = b 2 /N (see Fig. 3, where we have 
reported the outcome of several simulations performed for different values of N and b). The 
existence of a specific function (5(A) is equivalent to the standard assumption made in the 
scaling theory @J34|, i.e. to assuming that the logarithmic derivative of the conductance, 

dlnG 

" s dtaiV < 27) 

is a function only of G. The representation of the numerical results in terms of rj and G is 
given in Fig. 4. 

From the inset of Fig. 3, it can be seen that in the metallic regime, i.e. for A 3> 1, 
G ~ aiA + a , that is the conductance has the expected "ohmic" dependence on the sample 
length (A ~ 1/A 7 "). In the i] representation, this is tantamount to saying that 77 ~ — 1 + ao/G 
with ao < 0. In the localized regime, A C 1, Fig. 3 reveals an exponential decrease of the 
conductance with the chain length, G ~ b Q exp(— 61/A), a behavior which corresponds to 
V = \n(G/b ). 

Let us now discuss some approximate theoretical expressions of G and compare them with 
our numerical results. First of all, we want to derive a sufficiently accurate expression for 
the effective Lyapunov exponents in order to determine the conductance from the Landauer 
formula in the limit of large b and A" but finite A. The extensive studies performed in 



Ref. [|12] have suggested that 



7,(6, jV)« 7< (6, 00) -^-j- b , (28) 

is an expression convincingly tested in the bulk of the spectrum. Moreover, it was shown 
that 

7,(6,00)^—^ + -^ (29) 
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where £l(x) represents the asymptotic Lyapunov spectrum, while u(x) is the leading correc- 
tion term. By substituting Eq. ( ^9|) in Eq. ( ^8|) , one obtains 



7((6 , w) ^ + ^)__a_ ( 30) 

As the leading contribution to the conductance is given by the "small" Lyapunov exponents 
(i.e. x — I); ^ is convenient to expand 0(%) around x — 1- By retaining the leading terms 
(recall that 0(1) = 0) one obtains 

M ^ n <l^ + n JAgl + ^l^ (31) 

where Oi is the absolute value of the first of derivative of O(x), while 2 is the second 
derivative. 

From the Landauer expression (|14D , one can see that the argument of the hyperbolic 
cosines is 

„ , rWr 20i 20 2/ x o a;(l)-Oi 2A , . 

a, ee 2 7i (6, iV)iV « + -^(j - I) 2 + - ^ (32) 

where j = b — i + 1 Accordingly, one can see that in the limit b, N —>■ oo but finite A, the 
second and the last term in the r.h.s. vanish and thus we neglect them. 

Therefore, we are left with a linear expression in j, with two coefficients to be determined, 
namely Oi and ou(l). From the analytic knowledge of the minimum Lyapunov exponent 
(j = 1), we have that = 3/2. Moreover, previous simulations for the Lyapunov spectra 
have clearly indicated that the slope is Oi = 1.5. As a consequence, we are left with 

aj = 3j/A , (33) 

which, upon substitution into Eq. ( |TS ) implies 

oo 2 

G = §i + cosh(3j/A)' (34) 

where, in view of the limit b — > oo, the sum is extended to infinity. 

In the metallic regime, the conductance expression can be written in a compact form by 
transforming the sum into an integral, 

G = 2A r~ dy 

3 J3/A 1 + cosh?/ 

where we have introduced the variable y = Sj/A. The solution of the integral finally yields, 

G = 3(exp(3/A) + l)- (36) 



In the limit A ^> 1 one gets 



;oo 

G = ai A + a = - 1 (37) 



11 



with a% = 2/3 and ao = —1 . The first term in (|37j ) coincides with the theoretical result 
according to which the leading term for the conductance in the metallic regime is the ratio of 
the localization length to the size N of the sample. One should stress that the finite size 
corrections are neglected in this term (see, e.g., ||). A more important correction is given 
by the term ao and is related to weak localization effects. The meaning of this quantum 
correction is that the back scattering is larger in comparison to the classical result and, 
therefore, the conductance is smaller than the contribution given by the first term. In the 
theory (see the review || and references therein), the value of ao is known to be ao = — 1/3 
It is, however, interesting to note that a simple estimate of this correction |3~7[ gives 



the value a = — 1 as in our expression fl3"T|). 

In the metallic regime, the best fit of our numerical data with the expression (^) gives the 
values a\ ~ 0.74 and ao ~ —0.41 . While the value of a± is close to the theoretical prediction, 
the correction term a is somehow different. A discussion about the correspondence between 
numerics and theory is given below. 

We now go to consider the opposite limit of strong localization, A 1. In this case, the 
sum ( |34] ) is dominated by the first contribution, 

G = b exp (-^j (38) 

with b = 2 and b\ = 3 . Let us first discuss the leading dependence given by the term 
6i . One should first stress that the coefficient b\ is different for different definitions of the 
average. Specifically, if one treat G as the average conductance (G), the value of &i was 
found to be b\ = 3/4. For the geometric average G , instead, (also known as the "typical 
conductance"), the theoretical prediction is b\ = 3 (see |38|,|6||). Therefore, the expression 
(|38|), which refers to the latter case, gives the correct exponential dependence. 

Less is known about the correction to the exponential dependence which in fl38|) is ab- 
sorbed in the term bo . To the best of our knowledge, this correction is unknown for G , the 
quantity of our numerical interest. At the same time, for the average conductance (G) the 
correction has the form bo ~ A 3 / 2 [3S| . One can suggest, however, that for the geometric 
conductance G the correction has the form b ~ A u with some value v > 1. Our data 
do not allow to estimate the correction to the exponential dependence. The best fit to the 
expression ( p8[) with &i as a constant, gives the following values, &i ~ 2.9 and b ~ 5.58. 
One can see that the value of b\ is quite close to the theoretical one, b± = 3. 

A comparison of the theoretical expression (|34|) with the data in a broad interval of the 
localization parameter A is presented in Fig. 3. Considering that no free parameters are 
present in the theoretical formula, the overall agreement represents a nice confirmation of 
the various approximations behind Eq. (|34]). Nonetheless a small but clear deviation can be 
seen in the intermediate region A ~ 1. It is, therefore, natural and important to understand 
whether the disagreement is to be ascribed to the misjudgement of some correction term, 
or to finiteness of b and N in the simulations. Moreover, it is important to recall that 
conductance fluctuations are not taken into account in the derivation of Eq. (|34]), since it is 
expressed in terms of the mean Lyapunov exponents. This analysis goes beyond the scope 
of the present paper. 

One can also describe our data for the conductance G by making use the phenomeno- 



logical expression suggested in [Id], 
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(39) 



with four parameters q. By fitting directly all the parameters, one finds an overall excellent 
correspondence with our data, but a careful inspection of the behavior in the metallic and 
localized regimes is less convincing. On the other hand, one can determine the parameters 
Cj, in such limit cases by exploiting the correspondence between q and a , a\ , b , b\ , see 
Eqs. (|37D and @, 

ci + c 2 c 2 c 4 C1 + C2 , . 

ai = ; a = — ; h = c 3 ; b = a (40) 

Unexpectedly, the resulting dependence (|39| ) with these (different) values of q, turns 
out to be practically indistinguishable from the analytical expression (|34]) . Therefore, the 
dependence (|39|) has no advantage in comparison with Eq. fl34]) . It should be noted that in 
||16|| , where the same problem has been studied with a different method (scattering matrix 
approach) and different boundary conditions (the coupling matrix elements have been chosen 
to be constant as in the leads, instead of random as in the bulk), the parameters q and 
the behavior in the metallic regime are slightly different. At the moment, we cannot give 
a convincing explanation of this discrepancy. In fact, theoretically, the influence of the 
coupling is still not understood, while, numerically, the quality of the data is not such to 
allow drawing a definite about the influence of the coupling. 

We have also performed a detailed analysis of the distribution V(G) in various regimes. 
In particular, it was found that the shape of the distribution is controlled by the same 
scaling parameter A. In the localized region, the distribution is approximately log-normal 
(see Fig. 5). The variance of InG decreases upon increasing A i.e. moving towards the 
metallic regime and in the localized regime satisfies the following relation 

var(\nG) = -A \nG - B 0l (41) 

where Aq ~ 1.98±0.03 and Bq ~ 3±0.4 (see Fig. 6). The value of the slope is fully consistent 
with theoretical expectations and with the numerical results obtained in ||16|| . However, let 



us recall that in the computation of Thouless conductance, it was found that Aq ~ 1 [15 



This means that fluctuations properties of Landauer and Thouless conductance are different 



(see also [1_7|). The value of the constant B is instead slightly different from what found in 
|T5| | . This should be again a consequence of the different connection between the disordered 
sample and the leads. 

Moving towards the metallic regime, the shape of the distribution V(G) becomes normal 
(see Fig. 7). This result is in perfect agreement with the theoretical predictions for quasi- ID 
conductors [pzj] . 

The last problem that we have addressed in this paper is related to the so-called Universal 
Conductance Fluctuations. As is known (see, for example, the review |J), the correct value 
of the variance of the dimensionless conductance is var(G) = 2/15. At the same time, 
direct application of the standard Random Matrix Theory to transmission matrices in the 
metallic regime gives var(G) = 1/8 ||. The difference is small but crucial, it indicates 
that the repulsion between the transmission eigenvalues is not correctly described by the 
logarithmic two-body interaction, which is the core of the Random Matrix Theory. In 
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this sense, it is important to study directly in numerical experiments how the fluctuations 
of the conductance depend on the degree of localization. We have performed simulations 
for different values of the scaling parameter A. Our data for the variance are reported 
in Fig. 8. As one can see, for A still approximately equal to 4.5, the variance remains 
close to 1/8, the value predicted by the RMT. However, a further increase of the scaling 
parameter up to A ~ 7.4 leads the fluctuations to approach 2/15, the value expected from 
either diagrammatic calculations or the diffusion-equation approach for transfer matrices 
(for the references see 0). One should stress that it is important to satisfy the condition 
6 <C A/" in order to distinguish the metallic regime from the ballistic one. The latter, highly 
non-generic, regime is reflected in the points in Fig. 8 which deviate from above when this 
condition is violated. The data clearly indicates that the variance reaches the limit value 
2/15 very slowly, an effect revealed also by a different approach, based on the scattering 
matrix determined through the numerical solution of the Lippman-Schwinger equation, see 



V. DISCUSSION 

In this paper we have studied the statistical properties of the Lyapunov spectra for 
quasi- ID disordered samples of finite size in the connection with electric conductance and 
its fluctuations. The key point for our study is represented by Landauer expression, giving 
the conductance in terms of the Lyapunov exponents of a suitable transfer matrix. The 
mathematical model we have used is the band matrix ensemble with random entries for the 
disordered part. 

The central question we have numerically investigated, is how statistical properties of 
the Lyapunov exponents are reflected in the properties of the conductance. The main theory 
which gives quite a good description of the properties of the conductance, is the standard 
Random Matrix Theory (RMT) applied to transfer matrices. The conjecture that the RMT 
is a good mathematical tool for a theoretical approach, is based on the observation that 
different channels are strongly but irregularly coupled. This approach has led to many in- 
teresting theoretical results which have found quite a good support in several numerical 
experiments. Nonetheless, some small deviations can be found with respect to the RMT 
predictions. The most famous concerns the value of the normalized variance of the conduc- 
tance fluctuations which is 2/15 instead of 1/8. The difference is small but important for 
the theoretical implications about the universal properties of the conductance fluctuations. 

The origin of the above discrepancy has been attributed to different structure of the 
Lyapunov spectra compared to the spectrum of truly random matrices. For instance, in 
the RMT, the eigenvalues occupy the infinite line (for infinite random matrices with a fixed 
variance of the matrix elements), while the Lyapunov spectrum of the product of transfer 
matrices is symmetric with respect to the origin, a circumstance yielding an additional (non- 
random) repulsion between negative and positive "levels", which is absent in completely 
random matrices. On the other hand, in the diffusive regime, given the large number of 
channels, such an effect may be neglected. Indeed, earlier numerical experiments have 
shown quite a good agreement with the RMT prediction; for example, the level spacing 
distribution was found to be well described by the famous Wigner-Dyson (WD) surmise (for 
the lowest Lyapunov exponent in the diffusive regime). 
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In our study, we have paid special attention to the whole Lyapunov spectrum, perform- 
ing some qualitative tests aimed at emphasizing deviations of the level statistics from the 
Wigner-Dyson form. By making use of the normalized variance of the level spacing distri- 
bution, we have found that the statistical properties of the lowest Lyapunov exponent (as 
well as of the second one) are, indeed, described very well by the WD form. However, we 
have also detected a clear deviation from the RMT predictions for all the other exponents. 
In particular, it was found that the repulsion is weaker than that predicted by the RMT. 
This fact is quite unexpected since one might conjecture that the strongest influence of the 
negative part of the spectrum occurs for the lowest Lyapunov exponent, not for the bulk. 
Indeed, the lowest Lyapunov exponent has a "random" repulsion from one side only; on the 
negative side, the repulsion is of a different "regular" nature, due to the symmetry. However, 
it is precisely the lowest Lyapunov exponent to exhibit the best agreement with the random 
matrix approach (see also f23|). Thus we are lead to conjecture that besides the symmetry, 



a further feature must be present which contributes to differentiate a transfer matrix from 
purely random one. Anyway, the deviations observed in the bulk of the Lyapunov spectrum 
represent a further more detailed evidence of the failure of the RMT which has to be added 
to the difference between 1/8 and 2/15 for the conductance fluctuation (see also discussion 
and references in 0). 

To study the level spacing distribution in the general case (when the distribution does 
not coincide with the WD form), we have used the phenomenological expression suggested 
in |28| (see details in [^]) for the description of the so-called intermediate statistics which 
occurs due to strong localization effects. As long as the effective repulsion ranges between 
(Poissonian) and (3 = 1 (WD-distribution for the Gaussian Orthogonal Ensemble, GOE), the 



expression could be compared with the Brody-distribution |3£| . However, in the present case, 
the effective repulsion is larger {(3 > 1), even diverging in the limit of strong localization, 
where the distribution becomes an increasingly narrow Gaussian (a feature already discussed 



in the literature P2"|). Our analytical expression for the level spacing distribution allows to 
extract the repulsion parameter and to determine its dependence on the degree of localization 
(see Fig. 2). In analogy to other applications ||3l|j2"8ll , the effective repulsion turns out to 



depend only on the scaling parameter b 2 /N. 

In the second part of the paper, we have investigated various properties of the conduc- 
tance in a wide range of the localization parameter A. In particular, we have numerically 
computed the dependence of the geometric ("typical") conductance G on A, from the very 
localized to the diffusive regime. One should recall theoretical results on the dependence of 
the conductance on A are available only in the limit cases of a strong localization and in the 
metallic regime. 

Starting from the Landauer expression and taking into account the main properties of 
the Lyapunov spectra |ll,[12|, we have obtained an analytical expression for the conductance 
with no free parameters. Evaluation of this expression in the two opposite (metallic and 
localized) limits have shown that this expression gives, in general, the correct behavior in 
comparison with some theoretical predictions. In particular, in the metallic regime the 
leading term is exactly the same as predicted by the theory, although the weak localization 
correction is slightly different. 

A clear (but not very strong) difference has been found in the intermediate region A « 1. 
The source of this discrepancy is still not clear, in any case, one should emphasize that 
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in this region, both the coupling to the leads and finite size effects can play a significant 
role (we would like also to stress that there is no analytical prediction in this region for 
G). Moreover, we would like to note that the question of the influence of the coupling 
on the statistical properties of the conductance is not theoretically understood even for 
the relatively simple case of quasi- ID geometry. Recent numerical experiments with 
the same model ([!]) have shown that by changing the degree of coupling, the structure of 
scattering states changes dramatically In our present study, we have used the so-called 
"matching" principle according to which the variance of coupling elements is the same as 
that one in the bulk (see also |HJ); in this case the influence of the coupling is expected to 
be negligible theoretical predictions can be used. 

Finally, our results for the conductance fluctuations are in excellent agreement with the 
theoretical predictions. In particular, the distribution of the conductance has been found to 
be log-normal in the localized regime and of Gaussian in the diffusive regime. Furthermore, 
in the localized regime, the proportionality coefficient between the variance of the logarithm 
of the conductance and the logarithm of the conductance itself, is very close to the theoretical 
value 2. It is interesting to note, that when determining the conductance according to the 
Thouless definition via the curvature of the levels of a closed systems, this value was instead 
found to be 1 [15]. This fact indicates that fluctuation properties of the curvature in the 



localized regime cannot be associated with that of conductance fluctuations (see also [|Tj 

In the metallic regime, we have carefully studied the universal conductance fluctuations. 
One should stress that, numerically, this problem is far from trivial, since the sample length 
N has to be larger than band size b (defining the mean free path) but smaller that b 2 
(N <C b 2 , where b 2 is a measure of the localization length in infinite samples). On the other 
hand, one also needs a sufficient statistics in order to discriminate between 1/8 and 2/15 for 
the value of the variance of the normalized conductance. We have performed an extensive 
numerical analysis to overcome the above mentioned computational difficulties. The main 
result we have obtained is the observation of a slow convergence to the theoretical value 
2/15 upon increasing the localization parameter A. Nonetheless, we have also found a "pre- 
diffusive" regime, A « 4.5 where the variance stays very close to 1/8. This is an indirect 
indication that many Lyapunov exponents come into play, besides the lowest one, which 
have fluctuations slightly different from that given by conventional RMT. Thus, both the 
universal conductance fluctuations and the origin of the difference between 1/8 and 2/15, 



predicted theoretically in [gj], are nicely confirmed by our numerical study of the quasi-lD 
model. 
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FIGURES 



FIG. 1. Parameter y (see Eq. (p3|)) for different values of the scaling parameter A = b 2 /N 
versus the normalized number \ = i/b °f t ne Lyapunov exponent. The horizontal line corresponds 
to the Wigner-Dyson surmise. 

FIG. 2. Effective repulsion parameter (5 is plotted versus A for the minimum Lyapunov expo- 
nent 7b . 

FIG. 3. Average logarithm of the conductance versus the logarithm of the scaling parameter 
A for different values of the parameters N, b with N = L + 26. The solid curve is the expression 
(31). In the inset the geometrical average of conductance versus A is given for the metallic regime 
(Zoo ~ b 2 3> N). The best least square fit gives y = ax + b with a ~ 0.75 and b ~ —0.6. 

FIG. 4. Scaling function r/. The horizontal line gives the asymptotic limit r] = 1 for G ^> 1. 
Another line is used as a guide for the asymptotic slope of rj in the localized regime. In the inset 
we present the data for r] in the G » 1 limit. Straight line corresponds to the best least square 
fit (see in the text). 

FIG. 5. Log-normal distribution of the conductance in the localized regime. 

FIG. 6. Variance of the logarithm of the conductance against the average logarithm of the 
conductance in the localized regime A <C 1. The fit shows that the variance of the logarithm is two 
times larger than the mean value of the logarithm of conductance. 

FIG. 7. Normal distribution of the conductance in the metallic regime A> 1. 

FIG. 8. Variance of the logarithm of conductance versus the logarithm of the scaling parameter 
A. The low straight line is 1/8, the upper is 2/15. The error bars are of the order of the symbols. 
For all cases, we have averaged over more than 2500 different realizations of the disorder. 
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